Where two fractals meet: the scaling of a self-avoiding walk on a percolation cluster 
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The scaling properties of self-avoiding walks on a d-dimensional diluted lattice at the percolation 
threshold are analyzed by a field-theoretical renormalization group approach. To this end we recon- 
sider the model of Y. Meir and A. B. Harris (Phys. Rev. Lett., 63, 2819 (1989)) and argue that via 
renormalization its multifractal properties are directly accessible. While the former first order per- 
turbation did not agree with the results of other methods, we find that the asymptotic behavior of a 
self-avoiding walk on the percolation cluster is governed by the exponent v p = l/2+£/42 + 110£ 2 /21 3 , 
e = 6 — d. This analytic result gives an accurate numeric description of the available MC and exact 
enumeration data in a wide range of dimensions 2 < d < 6. 

PACS numbers: 64.60.Ak, 61.25.Hq, 64.60.Fr, 05.50.+q 



Polymers and percolation clusters are among the most 
frequently encountered examples of fractals in condensed 
matter physics 0, 0, ■ When a long polymer chain is 
immersed in a good solvent its mean-square end-to-end 
distance R 2 scales with the monomer number N as: 
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with the exponent vsAw(d) which depends on the (Eu- 
clidean) space dimension d only. This scaling of poly- 
mers Q is perfectly described by the self-avoiding walk 
(SAW) on a regular d-dimensional lattice and the 
fractal dimension of a polymer chain readily follows: 
^SAW = 1/fSAW- For space dimensions d above the upper 
critical dimension d up = 4 the scaling exponent becomes 
trivial: vsAw(d > 4) = 1/2, whereas for d < d up the 
non-trivial dependence on d is described e.g. by the phe- 
nomenological Flory formula [l| j^saw = 3/(d + 2). This 
found its further support by the renormalization group 
(RG) e — 4 — d-expansion known currently to the high 
orders Q: v SAW = 1/2 + e/16 + 15e 2 /512 + . . .. 

When a SAW resides on a disordered (quenched di- 
luted) lattice - such a situation might be experimentally 
realized studying a polymer solution in a porous medium, 
but is of its own interest as well - the asymptotic scaling 
behavior is a more subtle matter 0, 0, ■ Numerous MC 
simulations l a p, fill IllL Il2l fl3l and exact enumeration 
studies JI1 [H O^aTEHTmliSl, which last since early 
80-ies pf, lead to the conclusion that there are the fol- 
lowing regimes for the scaling of a SAW on a disordered 
lattice: (i) weak disorder, when the concentration p of 
bonds allowed for the random walker is higher than the 
percolation concentration ppc and (ii) strong disorder, 
directly at p = ppc- By further diluting the lattice to 
P < Ppc no macroscopically connected cluster, "perco- 
lation cluster" , remains and the lattice becomes discon- 
nected. In regime (i) the scaling law JIJ is valid with 
the same exponent vsaw for the diluted lattice indepen- 



dent of p, whereas in case (ii) the scaling law (QJ holds 
with a new exponent v p ^ ^saw- A hint to the phys- 
ical understanding of these phenomena is given by the 
fact that weak disorder does not change the dimension 
of a lattice visited by a random walker, whereas the per- 
colation cluster itself is a fractal with fractal dimension 
dependent on d: dpc(d) — d — /3pc/z/pc,where /3pc and 
v-pQ are familiar percolation exponents [2|. In this way, 
v SPw(d) must change along with the dimension dpc of 
the (fractal) lattice on which the walk resides. A modi- 
fied Flory formula Q for the exponent of a SAW on the 
percolation cluster i/ p = 3/ (rfpr: + 2) al ong with results of 
similar theoretical studies |2lj, I22L \'2'A |2J, |2jj|, l26t 1271 |2j| 
gives numbers in an astonishing agreement with the data 
observed (see Table QJ. Since d™ = 6 for percolation 0], 
the exponent v p (d > 6) = 1/2 |32| . 

Although the Flory-like theories [H |H |H |H 
HE H^i H3 offer good approximations for v p {d) in a 
wide range of d, even more astonishing is the fact that 
up to now there do not exist any satisfactory theoreti- 
cal estimates for v p (d) based on a more refined theory, 
which takes into account non-Markovian properties of 
the SAW, a task which was completed for regular lat- 
tices already in mid-70-ies 0. Existing real-space RG 
studies 0> EH EE HH give satisfactory estimates for 
d = 2, whereas the field-theoretical approaches aimed 
to describe the situation at higher dimensions lead to 
contradictory conclusions. In particular, the field theory 
developed in Ref. 01 supported d up — 6 and presented 
a calculation of v p in the first order of e = 6 — d. How- 
ever the numerical estimates obtained from this result 
are in poor agreement with numbers observed by other 
means, leading in particular to the surprising estimate 
i/p ~ ^saw m d = 3 (see Tabic [|J. In turn, a subse- 
quent study [iiil ] even questioned the renormalizability 
of this field theory and suggested another theory with 
d up = 4 which is obviously disproved by computer sim- 
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0.76(8) 
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0.595 


0.571 
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TABLE I: The exponent v p for a SAW on a percolation clus- 
ter. FL: Flory-like theories, EE: exact enumerations, RS, RG: 
real-space and field-theoretic RG. The first line shows vsaw 
for SAW on the regular lattice (d = 2 ^3|, d = 3 H3|)- 



illations and exact enumerations at dimensions d = 4, 5 

mm. 

There is another important reason, why the scaling of a 
SAW on a percolation cluster calls for further theoretical 
study. As it became clear now, higher-order correlations 
of a fractal object at another fractal lead to multifractal- 
ity |35|. Recently studied examples of multifractal phe- 
nomena are found in such different fields as diffusion in 
the vicinity of an absorbing polymer l3B , random resis- 
tor networks |37j |. quantum gravity |38| . A SAW on a 
percolation cluster is a good candidate to possess multi- 
fractal behavior. Indeed such behavior is found in com- 
puter simulations |l9|. moreover it naturally emerges in 
the RG scheme, as we will explain below. 

Let us consider a diluted lattice with sites x^ in terms 
of variables Pij = 0, 1 that indicate whether a given bond 
between the sites x^ and Xj is present or not. To describe 
the critical properties of SAWs on this lattice following 
the idea of de Gennes Q we introduce m-component spin 
variables S* Q (xi), a = 1, . . . , m, and evaluate the theory 



for m = 0. To allow for the averaging over the quenched 
disorder the spins are n-fold replicated which gives for 
the Hamiltonian: 

T ^ m n 

e-™» =< cxp{- - Y^Pij E E ^(x,)Sf (x,} > p 

i.j a = l /3=1 

(2) 

where we denote by < . . . > p the average over the random 
variables p^ which take the value 1 and with probabili- 
ties p and (1 — p) respectively, and K is an interaction pa- 
rameter. In the following we will work with a field theo- 
retical representation of the effective Hamiltonian defined 
in J2J. This is achieved [LJ| via a Stratonovich-Hubbard 
transformation to tensor fields ipk(x) with components 
^fc-/3i a /3 fc ( x ) conjugated to the product II* : =1 S'fj(x) of 
k components of the replicated spin with fj\ < ... < {3k ■ 
This results in the effective Hamiltonian up to order ip 3 

0: 

= l -l d d q^(r k + q 2 )Mn) ■ M~<l) + 
J k 

| J d d ^ 3 (x), (3) 

where ip^ (q) is the Fourier transform of ip^ (x) , the inner 
product reads: 

and ?/> 3 (x) is a symbolic notation for a product of three 
fields. Only those cubic terms ip 3 are allowed for which all 
pairs (cti, (3i) appear exactly twice. A second condition on 
the diagrammatic contributions to perturbation theory 
can be derived from the de Gennes limit m — 0, namely, if 
any index (a, /?) appears only on the internal propagator 
of a diagram, then its contribution vanishes. 

We note the unusual dependence of "masses" on k. 
This is reminiscent of the fact that in the m = limit 
the theory @ becomes multicritical (l4l [3!| . This has 
impact on the renormalization of the theory J3Jl as we 
will show in the following. 

We choose to calculate the critical properties of the 
theory by analyzing its vertex functions, in particular 
rW(q), T^({q}), and T^ 1 ) ({<?}) where the latter in- 
cludes an insertion of the ipitp operator. Each of these 
T- functions will depend on the family of masses {rk}- 
The Feynman graphs of the contributions to the two- 
point vertex function (q) in the two lowest orders are 
shown in Fig. ^ The contributions to r^ 2,1 ) are found 
from this by placing an insertion on each of the inner 
propagator lines. These integrals are evaluated then in 
dimensional regularization in dimension d = 6 — e and 
minimal subtraction |4o| using a Laurent-expansion in 
e. Usually the renormalization of the vertex functions 
is defined in terms of ^-factors in such a way that the 
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FIG. 1: The Feynman graphs of the vertex function V^ 2 \q) 
in the two lowest orders. 

products Z^ 2 \ Z W T {3 \ Z^TW are free of e-poles. 
However, the insertion of the ^:"0-operator together with 
the fc-dependence of the masses leads to the following 
renormalization procedure. The vertex function I^ 2 ' 1 ) 
even when evaluated at zero mass remains ^-dependent: 

r (2,l) = r (2,l) + fcr (2,l) + fc 2 r (2,l) + _ _ _ (4) 

and it can not be renormalized by one multiplicative Z- 
factor. The essential feature of this expansion is that each 
term shows a different scaling behavior. In this way the 
multicriticality r ecog nized already by Derrida [13 an d 
Meir and Harris |l4j manifests itself in our present for- 
malism and leads to a spectrum of exponents. Instead of 
a single Z-factor Z^2 a whole family of factors Z^.', is 

necessary to renormalize each T^ 2 ' 1 ^ in Q. This allows 
to define RG-functions 

= A in z (5) 

an 

that describes the RG-flow of the coupling with respect 
to the rescaling parameter k and 

VH = ^ In Z^ = A l n Z ^ (6) 

which define the anomalous dimensions of the corre- 
sponding operators. At the stable fixed point w* with 
/3(w*) = the family of correlation exponents is given 
by 

= [2 -»?(«;*) +r?^K)]- 1 . (7) 

We note that = vpc and v ^ = a s introduced 
above, whereas the i/W for i > 2 are connected with 
higher order correlations. The /3-function 10 is the fa- 
miliar RG function of the ip 3 Potts model |41| . 

The explicit calculations proceed as follows: (i) One 

(2) 

starts with the vertex function r^ fc corresponding to the 
propagator of the field ipk- (ii) For the masses one in- 
serts the expansion = /x^°^ Ujk° . (iii) The insertion 

fi (2) 

of is defined by the derivative -§j^\p k evaluated at 
zero mass for // = 0. (iv) Performing the summation 
over the replica indices the contributions to the different 

(2 1) 

' are generated by rearranging the expansion in k. 

(2 1) 

One finds the multiplicative renormalization for T t ' for 
appropriate linear combinations of the different orders of 
k. 



Following this procedure we obtain e = 6 — d ex- 

(i) 

pansions for rj) h . qh . Substituting them together with the 
known result HJ V = -e/21 - 206e 2 /21 3 into we ar- 
rive at the following spectrum of correlation exponents: 

z/ 0) = vpc = 1/2 + 5 e/84+ 589 e 2 /42 3 , (8) 
= is p = l/2 + e/42 + 110e 2 /21 3 , (9) 
v {2) = 1/2 + e/24+ 13907 e 2 /1100736, . . . . (10) 

By (JHJ we recover the familiar e-expansion for the perco- 
lation exponent i^p [4ll| and in © we extend the first 
order result for v v [l4|. The physical interpretation and 
properties of the remaining exponents v^ 1 ' of the fam- 
ily is the subject of a separate study 01- Contrary to 
the family of ^-exponents defined in Ref. 0], the i/W 
govern the non-trivial scaling of properly defined cumu- 
lants of the distribution of SAWs for given end-to-end 
distance 01 • Evaluating the result for is p (0 by direct 
substitution of e — 6 — d one finds nearly perfect cor- 
respondence with available MC and exact enumeration 
results over the range d = 2, ... ,5, see Table H] This 
presents a qualitative improvement over the linear result 
as seen in Fig. [21 where we also show that the result is 
in between the limits given by the shortest and longest 
SAWs on percolation cluster |44j . 




FIG. 2: The correlation exponent v p . Bold line: Q, thin line: 
one- loop result filled boxes: Flory result u p = 3/(dpc + 2) 
with dpc from |43|. Exponents for the shortest and longest 
SAW on percolation cluster |44| are shown by dotted lines. 

A rather peculiar finding is that results of the phe- 
nomenological Flory-like formulae evaluated using the 
fractal characteristics of the percolation cluster are nu- 
merically very close to our result in the same region of 
dimensions. Note however the ambiguity fljj in defin- 
ing a Flory-like scheme leading to the different results in 
Tabled 

The ip 3 theory as applied to the present problem in- 
evitably has the upper critical dimension d np = 6. This in 
particular allows us to describe the discussed non-trivial 
scaling for dimensions d = 4, 5. This is out of reach 
following the approach of Ref. 34] which gives trivial 
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scaling for d > d up — 4 and relies on a </> 4 -theory with 
two couplings of different symmetry. Moreover, in the de 
Gennes limit m = the symmetry of the two couplings 
coincides leading back to a theory of SAWs on the pure 
lattice with a redefined coupling parameter, a fact nei- 
ther exploited in Ref. [34| nor in the similar approaches 

mm. 

From the physical point of view, our result for the ex- 
ponent 2/ p together with the data of EE and Flory-like 
theories (see Table [|J predicts a swelling of a polymer 
coil on the percolation cluster with respect to the pure 
lattice: v v > i^saw for d = 2 — 5. Up to now, this 
phenomenon has clearly been observed only in MC sim- 
ulations for d = 2 12]. Although simulations on d = 3 
percolation clusters have been claimed to show this effect 
Hi El El i these studies were subsequently criticized 
for using inappropriate data analysis 0, E3, El and for 
lack of accuracy. At d = 3 our formula (jlOJI predicts a 
13% increase of v v with respect to ^saw which is larger 
than at d = 2 (5%) and should be more easily observed by 
current state-of-art simulations. Given that even at d = 2 
we are in nice agreement with MC and EE data and the 
reliability of the perturbative RG results increases with 
d, this number calls for verification in MC experiments 
of similar accuracy. 
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